Coulomb Effects on Electron Transport and Shot Noise in Hybrid 
Normal-Superconducting Metallic Structures 
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We analyze the effect of electron-electron interactions on Andreev current and shot noise in 
diffusive hybrid structures composed of a normal metal attached to a superconductor via a weakly 
transmitting interface. We demonstrate that at voltages/temperatures below the Thouless energy 
of a normal metal Coulomb interaction yields a reduction of both Andreev current and its noise by 
a constant factor which essentially depends on the system dimensionality. For quasi-ld structures 
this factor is ^ ' where A^ch and g are respectively the number of conducting channels and 
dimensionless conductance of a normal metal. At voltages above the Thouless energy the interaction 
correction to Andreev current and shot noise acquires an additional voltage dependence which turns 
out to be a power-law in the quasi-ld limit. 
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I. INTRODUCTION 

It is well established that low temperature charge 
transport through an interface between a normal metal 
and a superconductor (NS) is dominated by Andreev 
reflection^: An electron with energy below the super- 
conducting gap A enters the superconductor from the 
normal metal, forms a Cooper pair together with another 
electron, while a hole goes back into the normal metal. In 
addition to Andreev reflection electrons can be scattered 
at the interface potentials and/or impurities. The combi- 
nation of "normal" and Andreev reflection mechanisms is 
essential for proper understanding of transport phenom- 
ena in proximity structures composed of superconduct- 
ing and normal metals. This applies both to equilibrium 
phenomena, such as dc Josephson effect in diffusive SNS 
hybridsSiiii^ and Meissner effect in NS systems^i^, and 
to non-equilibrium effects, such as dissipative transport 
of subgap electrons across NS interfaces^. 

For weakly transmitting NS interfaces the correspond- 
ing Andreev conductance Ga turns out to be rather small 
as a second order effect in the interface transmission. At 
the same time interplay between disorder and interfer- 
ence effects in the normal metal may strongly enhance 
Andreev conductance at sufSciently low energiesiSiiiiiSii^ 
leading to the so-called zero-bias anomaly (ZBA) on the 
current-voltage characteristics. In the case of diffusive 
metals Andreev conductance grows as 



Z{uj) = R Coulomb interaction yields a power-law ZBA 
which, being combined with Eq. Q, yields^'* 



Ga(xI/Vt, Ga Oil/ W 



(1) 



with decreasing temperature and applied voltage V. 

What is the effect of Coulomb interaction on the An- 
dreev conductance of NS structures? In the limit of 
weakly transmitting NS interfaces this question was ad- 
dressed in Ref. within a simple capacitive model and 
the Coulomb blockade of Andreev reflection was pre- 
dicted at sufSciently low energies. Huck et al^ stud- 
ied the problem by modeling the effect of interactions 
by means of an effective electromagnetic environment 
with an impedance Z[Ltj). Similarly to the the case of 
normal tunnel junctions^^' , for the Ohmic environment 



Ga oc t8/9-i/2^ Ga oc v'/^-^/^ 



(2) 



respectively for eV < T and eV > T. Here and below 
we define the dimensionless conductance as g = Rq/R, 
where Rq = h/e^ is the quantum resistance unit. 

While the approach''^^ takes into account Coulomb in- 
teraction within a phenomenological scheme of the so- 
called P(£^)-theory^^-^^, the question remains if and un- 
der which conditions this scheme is sufficient to account 
for the effect of electron-electron interactions in a dis- 
ordered normal metal attached to a superconductor. In 
this paper we develop a microscopic analysis of this prob- 
lem and evaluate the interaction correction to Andreev 
conductance of diffusive NS structures. We demonstrate 
that both the form and the magnitude of this interac- 
tion correction essentially depend on the dimensionality 
of the sample as well as on the relation between tempera- 
ture and applied voltage to the relevant Thouless energy 
of the structure. In particular, we predict that at suffi- 
ciently low energies the interaction correction saturates 
to a constant value, thus providing a temperature- and 
voltage-independent renormalization of Andreev conduc- 
tance. This saturation effect was recently observed^^ in 
NS structures composed of a superconductor and multi- 
walled carbon nanotubes (MWNT). 

The structure of the paper is as follows. In Sec. II 
we will deflne the model for the NS system and outline 
the formalism of a non-linear a-model which will be used 
for our analysis. Sec. Ill contains the derivation of the 
general expression for the part of the action responsible 
for Andreev processes in our system. With the aid of 
this expression in Sec. IV we will flnd the Andreev cur- 
rent in the presence of electron-electron interactions. Its 
analysis in a number of important limiting cases is pre- 
sented in Sec. V. In Sec. VI we will establish a general 
relation between the current and shot noise in the pres- 
ence of electron-electron interactions. In Sec. VII we 
will briefly discuss our key observations and compare our 
results with recent experimental findings—. 
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© is taken over the space and time variables and it is 
accompanied by the summation over matrix indices. 

The term Sr Q follows directly from the Kupriyanov- 
Lukichev boundary conditions^ which account for tun- 
neling of electrons between N- and S-mctals. Spatial in- 
tegration in this term is restricted to the insulating in- 
terface, subscripts =F denote matrices taken at the left- 
and right-hand sides of the interface. Other matrices in 
Eq. © are defined as follows 



FIG. 1: NS hybrid structure under consideration. 



II. MODEL AND BASIC FORMALISM 

Throughout the paper we shall consider the hybrid 
structure depicted in Fig. 1. A piece of a normal metal 
of a simple rectangular shape and of length L is attached 
to a bulk superconductor via an insulating interface and 
to a big metallic reservoir N' via a highly conducting in- 
terface. The transversal dimensions of the normal metal 
are Ly and Lz, hence, the cross-section of the NS inter- 
face r is just the product F = LyL^. In what follows it 
will be convenient for us to assume that the transmission 
of a tunnel barrier at this interface is sufficiently low, 
so that its resistance Rt exceeds that of a normal metal 
R = L/aV, 



Rt > i?. 



(3) 



where a = 2e^DNo is the standard Drude conductivity 
with e being the electron charge, D = vpl/3 standing for 
the diffusion coefficient, and iVo representing the density 
of states per spin direction at the Fermi surface. In ad- 
dition we will assume that the normal metal is shorter 
than the dephasing length L^, 



(4) 



Our theoretical analysis is based on the Keldysh rep- 
resentation of a non-linear cr-modeliSiiS which is formu- 
lated in terms of the effective action 



f , 
f , 



1 

-1 



(9) 



A{R, t) 



A 
A 



A = 



A(i2,t) 

-A*{R,t) 



The operator d acts as 
c 

Vg{t, t') - i-A{t)fzg{t, t') + g{t, A(t')f,.(10) 

c c 

Here t) is the vector potential. In Eq. (|10|) spatial 

arguments are suppressed for brevity. 

Finally, the V term depends on the Hubbard- 
Stratonovich fields Vi and V2 defined on the two branches 
of the Keldysh contour. These fluctuating flelds emerge 
after the standard decoupling procedure in the term de- 
scribing electron-electron interactions. We define 



V{R,t) 



v+l \v-i 
W^i v+i 



(11) 



where i is 2 X 2 unity matrix and — (V^i+V2)/2, — 
Vi — V2 define respectively "classical" and "quantum" 
fiuctuating fields. If an external voltage V is applied to 
the system, it should simply be added to the classical 
component V'^ . 

It is necessary to supplement the action de- 
scribing fermionic degrees of freedom, by the action for 
the bosonic fields . The latter action determines the 
correlators for these fields which will be specified below. 



S — ^0 + Sy, 



(5) 



where 



Sq = — — Tr 



St 



D {dgf + Ai {ifA + A - ey) 5 (6) 



Ae'^RtV 



^^vg-g^ 



(7) 



Here g{R, ^1,^2) is the 4x4 matrix that depends on one 
coordinate and two time variables. It obeys the normal- 
ization condition 



j dt'g{R,h,t')g{R,t' ,t2)^5{ti-t2). 



(8) 



The product of matrices g in Eq. © should be under- 
stood as a convolution, cf. Eq. ©. The trace in Eq. 



III. ANDREEV ACTION 

Let us first define the matrices g± in the ab- 
sence of both electron-electron interactions and tunnel- 
ing between metals. For the normal metal one has 
g+{R,t,t') = go{t,t'), where 



'"-^ g^ 



go{t,t') 



go 

2iT 



sinh[7rT(t - t')] 



(12) 



The Fourier transform of the Keldysh matrix 5^ reads 
g^{e) = 2 tanh[e/2r]f2. In what follows we will consider 
the values of voltages and temperatures much smaller 
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than the superconducting gap eV,T <C A. Hence, it will 
be sufficient for our purposes to specify the matrix g_ 
for a superconducting electrode only at energies much 
smaller than A. Under this condition we have 



-i 

1 


i 



S{t-t'). 



(13) 



Let us turn on electron-electron interactions. While in 
a large superconducting electrode their influence on the 
matrix g- can be safely neglected, the matrix g+ in the 
normal metal gets modified asi^ii^ 



g+{R,ti,t2) = e 



5o(ti,t2)e^^(«'*=). (14) 



Here the matrix K has the structure 



K = 



(15) 



and the phase fields ip^ satisfy the equations 

{Dq^ + iuj)ip- {q,Lo) + eV^ {q,uj) = 0, 
[Dq^ -iuo)^+{q,u)-eV+{q,u) = 

~Dq^Lp-{q,uj)coih^. 



(16) 



One can easily check that the part of the variation of 5o 
^ linear both in ^-variations and in the fields van- 
ishes under the transformation l|14|l - l|16|l . These equa- 
tions are sufficient to account for the effect of Coulomb 
interactions for the problem in question. 

Now let us include tunneling between metals into con- 
sideration. In order to eliminate the dependence of the 
matrix g+ on the fluctuating fields let us permute the fac- 
tors exp(±iif) in the term St (0 to act on the matrix 
5_ . After this transformation we obtain 



;Trr Q-.9+, 



(17) 



where 



B{R,t) A{R,t) 



A = 



B = 



ie '^^'^ cos 



(18) 



sm ip 







What remains is to evaluate the deviations of g^ from 
its normal state value H12|) due to tunneling of Cooper 
pairs from the superconducting electrode into the normal 
metal. It is convenient to employ the following parame- 
terization 



Q — iig+u, u — u ^ — 







iT 



sinh[7rT(t-t')] 



(19) 



which brings the matrix Q to the diagonal form in the 
absence of the proximity effect. The latter effect can 
can then be interpreted in terms of fluctuations of Q 
parameterized as 



Q 



with 



(1 + W + WV2 + ...) 



1 

-i 



(20) 



(21) 



The matrix W is in turn parameterized via the Pauli 



matrices Tx 



(22) 



Wlfj; + WiTy WqI -I- WsTz 
Wo i -t- W-iT^ WiTx + W2Ty 

The functions Wi , and Wi can be expanded in the cigcn- 
functions ijj of the diffusion equation 



Wi{x,y,z) ^w.t{q)tpq(x,y,z), 



(23) 



obeying the boundary conditions dip/dxlx^o = Oi 
'0|x=L = and impenetrable boundary conditions at 
y = 0,Ly and z — 0,Lz. We choose 



ipq{x,y, z) oc cos qxX cos qyy cos q^z, 
I' 



l,m,n = 0, 1, 



(24) 



The normalization constant is determined by the condi- 
tion / dxdydz'ipq{r) = 1. 

Using the effective action ^ , we can find all non-zero 
averages for the products of the functions Wi, and Wi. 
Below we will only need the following averages^^ 

(wi(q, ei, e2)?i;i(-q, 63, 64)) = 

1 (27r)^(5(ei - e4)(5(£2 - £3) 

ttNo Dq^ - i{ei + £2) 
{Wi{q,€i,€2)wi{-q,e3,e4)) = 

(27r)25(ei - e4)<5(e2 - ea) 



1 



Dq^ +i{ei +62) 



(25) 



which correspond to the Cooperons (i — 1,2). 

Now we are ready to evaluate the contribution of An- 
dreev processes to the effective action. In the weak tun- 
neling limit it is sufficient to expand exp(iS') up to the 
second order in St and then to average the resulting 
expression over the fermionic degrees of freedom with 
subsequent re-exponentiation of the result. In the limit 
Rt S> R^ the amplitude of the anomalous Green function 
penetrating into the normal metal remains much smaller 
than unity. Hence, in Eq. H2U|I it is sufficient to retain 
only the linear terms in W. The resulting Andreev con- 
tribution to the action then takes the form 

(Trr,C-(ri)*.f.>V(r,)-Trr,C-(rg)*.f.H'(r,)). 
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The space integration is performed here along the same 
interface F, but independently for the terms containing 
Fi and F2. The matrix Q_ is defined by the relation 
Q_ = uQ^ii. Employing Eqs. H18|) . Ijiyi) and performing 



the averaging with the aid of Eqs. H25|l . we arrive at 
the final expression for (5S'^[(/3*] which will be extensively 
used in our subsequent analysis. We find 



Dq^ 



tanh - tanh ^ ' ' 



2T 2T 
+ sin [2ip+{r,t) - 2ip+{r',t')] cos[v5" (r, i)] sin[.^~ (r', <')] 
2iE 



coth ile-2^(^"'('''*)-'^^('''^*'» sin[v7- (r, t)] sin[(^- (r', t')] 



(Dq^f+AE"^ 

Dq^ 
{Dq^f+AE^ 



, E+{uj/2) , E-{uj/2) 

tanh ^—LJ. + tanh ^ ' ' 



2T 2T 

e2'(^^ cos[</j-(r, t) - (p-(r',i')] 



in [2ip+{r,t) - 2ip+{r',t')\ cos[(^"(r,t)] sin[(p~(r',t')] 



(27) 



IV. ANDREEV CURRENT IN THE PRESENCE 
OF INTERACTIONS 

Let is first evaluate Andreev current across the NS in- 
terface. Substituting Eq. (|27|l into the formal expression 
for the current 

I 



and performing averaging over the fluctuating fields ip- 
we arrive at the general result for the Andreev current 



^ _ ttRD 



2eRlTL 
Dq^ 



dEduj 



tanh 



E + eV -{uj/2) 
2T 



tanh 



E~eV +{uj/2y 
2T 



[Dq^ 



iE^ 



P{io,r,r') 



1 - e-"^^/^ 

I _ g{uj-2eV)/T- 



(29) 



Here the function 



^J{t,r,r') ^ /2i<p2{r,t)~2iipi{r'fi) 



(30) 



accounts for the effect of interactions described by the 
fluctuating phase fields ipi^2 ~ 'P^ ± (<<5 /2) defined on 
the two branches of the Keldysh contour. Averaging in 
Eq. H3U|) is performed in the absence of the external volt- 
age. Note that the form of this function is reminiscent to 
that used in the standard P(£')-theoryi^4i& and the result 
for the Andreev current (|29|l looks somewhat similar to 
the expression derived within the framework of the latter 
theory in Ref. |3 In contrast to the latter work, how- 



ever, our approach and resulting Eqs. (^51) . l|Hn)l fully 
account (i) for electron coherence across the whole NS 
structure and (ii) for the spatial dependence, introduced 
both by the functions ipq{r), describing fermionic fluctu- 
ations, and by the function Jit, r, r'), describing bosonic 
(phase) fluctuations. 



These fluctuations depend on the dimensionality of 
the normal sample. For our purposes it will be suf- 
ficient to employ the standard random phase approx- 
imation (RPA) and to express the correlators of the 
fiuctuating Hubbard-Stratonovich fields via an effective 
d-dimensional dielectric function e(w, k) of our system. 
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One finds2a 

{eV+{t,r)eV+{0,0)) 



{eV+(t,r)eV-{0,0)) = 



-Im(y*(fc,c.))coth(^)e— 



dtud'^k ^ 



— iLot-\-ikr 



(27r)« 

(ey"(t,r)ey"(0,0)) = 0. 
Here we have defined 
V*{k,u;) 



Vo-'{k) + ^[e{k,c,)~iy 



(31) 



(32) 



where the term Vo{k) stands for the unscreened Coulomb 
interaction: 



Vo{k) = 
Vo{k) - 



Vo(fc) =e^ln 1 + 



27re^ 
1 

ffc2 



d = 3, 
d= 2, 
<i= 1. 



(33) 



In the latter equation we have assumed Ly ^ 

^/T. In the three-dimensional case we obviously have 

V*ik,Lo)^Vo{k)/e{k,uj). 

Note, that in the case of lower dimensions d = 1,2 the 
term Vf^^{k) also accounts for the electric field energy 
outside the conductors. The corresponding interaction 
can be screened, provided there are metallic gates near 
the film/wire. For instance, in the case d = 2 one has 



Vo-{k) = ^ (1 



(34) 



where b is the distance between the film and the gate 
electrode. 

In order to evaluate the average for the phase fields 
(|30|l one can combine Eqs. (|16() . H31() with the standard 
Drude form of the dielectric function for the normal metal 



e(fc,w) = 1 + 



-iuj + Dk'^ 



(35) 



The resulting expressions for the function H30|l ob- 
tained in this way remain applicable at frequencies above 
the corresponding Thouless energy of the sample but 
should fail at lower frequencies, in which case discrete 
wave vectors ^ 1/L become significant. In order to cure 
this complication it is useful to depart from the repre- 
sentation of plane waves in Eq. (|31|l and to employ a 
different set of the (volume-normalized) basis functions 
^m{f) obeying the equation 



(36) 



with the corresponding eigenvalues Am. At the impene- 
trable interfaces the functions flmi'i') satisfy the bound- 
ary condition V„f2m(T") = 0. One should also require the 
functions Qrn{f) to vanish at the interface with a large 
metallic reservoir x = L. 

The corresponding generalization of the dielectric func- 
tion H35|) can be achieved by evaluating the response 
of a finite sample, characterized by the current density 
ji(r) = (jE{r) — D\7p{r) with E and p standing respec- 
tively for the electric field and for the charge density. One 



getsSi 



£af3{r, r', uj) = Saf3S{r - r') + 



(37) 



<yD {So^fidir - r') - VaVpd{r, r' , uj)) , 



where 



d(r,r',w) = ^ 



^r,i{r)^7n{r') 



(38) 



Implementing these modifications also into Eqs. H31() 
and employing Eqs. Ijlbfl we arrive at the following gen- 
eral result 



{DXl 



J(t,r,r')-45^ I ^Im 

(coth^coswt — isinutj VL„i{r)VLm{r') — (39) 
^coth^ (an(r)f^m(r) + Q.„,{r')a,n{r')) 



The latter formula results in J{t ~ {i/T),r,r') = 
J{—t,r,r') leading to the relation P{—E,r,r') — 
e~^/'^ P[ijj,r,r') as in the standard P{E) theory. Ac- 
cordingly, at T Eq. reduces to 



DR 



dl 



dV 2ttRITL 



2\eV\ 



^ j dE j d^rd^r'ijq{r)ijjq{r') 



PiE,r,r') 



Dq^ 



{Dq^Y + {2\eV\ - Ey 



(40) 



At frequencies exceeding the sample Thouless energy it 
is possible to perform averaging of Eq. H39|l over fluctu- 
ations depending on r -I- r'. Then one gets 



Jit,\r\) 



dwd'^k 
(27r)"+i 



Im 



V*{k.u)) 



coth — (cos ix}t cos fcr — 1) 



(L>fc2 ~ iujf 
- i sin Lut cos kr 



.(41) 



V. RESULTS 

A. Non-interacting limit 

Let us first reproduce the results for Andreev current in 
the non-interacting limit. In this case we have P{uj, r) = 



6 



2tt6{uj) and the standard expressions for the I — V curve 
are readily recovered. In the Hmit eV ^ T one finds 



dl 
dV 
dl 
dV 



R_ 

R 
2if 



if \eV\ < exh; 



\evy 



if \eV\ > exh. 



(42) 



Here and below exh — D/L? is the Thouless energy. 

In the linear regime eV <^ T for the Andreev conduc- 
tance Ga = dl/dV we obtain 



^a'' - if r < exh, 



(43) 



G 



(0) 



(23/2 _ 1)^(3/2) R 



40F 



RlV T 



if T » exh, 



where (2^/2 - l)C(3/2)/4v^ « 0.67. 

Now let us turn to the effect of electron-electron in- 
teractions. Both the value and the form of the interac- 
tion correction to the Andreev conductance essentially 
depend on the effective dimensionality of fluctuations of 
the phase fields. Depending on the voltage and temper- 
ature the leading contribution to the interaction correc- 
tion is given either by the modes with l,m,n ^ 1 (cf. 
Eq. (|24ll ) or by the modes with m (or n, or both) being 
equal to zero. Correspondingly, below we will distinguish 
between effectively 3d-, quasi-2d- and quasi-ld-type of 
behavior of the interaction correction. 



B. Effect of interactions: Bulk limit 



We start from the case of 3d fluctuations. For conve- 
nience we shall assume L ^ Ly ^ and consider the 
limit eV > T. At sufficiently high voltages eV > D/Ll 
from Eq. 1411) (with d = 3) we obtain 



J{t,r) 



cos{2Dtuj/r'^ 



e ^ cos VoJ ) div 



One observes that the relative "weight" of the sec- 
ond term in the brackets of (|47|) grows with voltage as 
{1/ g)y^eV/e^\i. This form is consistent with the stan- 
dard square-root energy dependence of the interaction 
correction to the density of states. In the second term 
in Eq. H47|l this dependence is exactly compensated by 
ZBA (|42f) resulting in the voltage-independent contribu- 
tion ^ R/ gR^. This compensation will not occur in other 
limits to be considered below. 

The value of k in the bulk limit can be estimated as 



l/r 



duj 

7^ 



{pFlf 



(48) 



Here t = l/vp is the elastic scattering time. 

While the result (|47|l holds at high voltages, at lower 
values eV ^ exh - as it was already discussed above - 
the Andreev conductance saturates to a constant value 



dl 
dV 



Re- 



(49) 



with the voltage-dependent correction to this formula of 
order {eVY / gei^^^. 

Note that in the bulk limit the net effect of electron- 
electron interactions is parametrically small and negative 
in the admissible range of voltages eV < 1/t. Indeed, 
having in mind that our calculation is performed in the 
limit ppl ^ 1 and making use of the estimate (|48|l . in 
Eq. H47|) one can expand exp(— k) to the first order in k 
and observe that the corresponding negative correction 
exceeds the positive one by a factor ~ 1/ y/\eV\T. 

Let us also point out that in order to obtain the An- 
dreev conductance in the voltage range eV <C T it is 
sufficient to simply substitute T instead of eV. As in 
the non- interacting limit (cf. Eqs. H42|l and (|43|l 'l. this 
procedure will yield the correct form of the interaction 
correction up to an unimportant numerical prefactor of 
order one. 



-erf 



sgn t — K, 



where erf(a;) = 2 cxp{—t^)dt/^ and 

V*{k,uj 



K = 4 



'0 

dujd'^k 
72^ 



Im 



{Dk^ - iLuy 



signcj. 



(44) 



(45) 



Making use of the condition P{u!,r) — for w < 0, we 
have for w > 

P(w,r) = e-''2T:d{Lu) + 
2e^e-'' 



TTUTUJ 

Thus we obtain 

dl Re 
dV 



_ g-r^W^cos (^ry/uj/2D^ . (46) 



2i?2 



3.64r 



\eV\ gL^ 



(47) 



C. Effect of interactions: Quasi-2d structures 

Let us now turn to low-dimcnsional fluctuations. We 
first consider systems which can be treated as quasi- 
2d ones, i.e. we assume that both the length and the 
width of the normal metal strongly exceed its thickness 
L ^ Ly ^ Lz- We shall address two different physical 
situations corresponding to the absence and to the pres- 
ence of massive metallic gate electrodes in the vicinity of 
the sample. 

In the absence of the gate electrodes we can use the 2d 
expression for Vo(fc) combined with Eq. We 

first evaluate the constant k which now consists of two 
contributions. 



K = K3 + K2, 



(50) 
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where the first one is again determined by high frequen- 
cies (w > D/Ll) and small scales (cf. Eq. if^ ). while 
the second one emerges from integrating over the param- 
eter region 27r crL^g ^ u and Dq^ <C w in Eq. for 
smaller frequencies uj ^ D/L^ 



K.2 = 



8Ly Ly loqT 

— - In — In , 

D ' 



ngL 



D 



(51) 



It is obvious from the above expressions that the con- 
tribution K3 exceeds K2 for thicker films, i.e. as long as 
Lz remains bigger than the mean free path I. In the op- 
posite case, however, the contribution of 3d fluctuations 
K3 is negligible, and in the expression for K2 one should 
substitute l/r instead of D /L^ in the arguments of the 
logarithms. Then one finds 



2L,, 



...-^Iniiln^. 
TrgL Dt D 



(52) 



We observe that the above expressions for the renor- 
malization constant k, contain double logarithmic fac- 
tors which are rather typical for 2d structures. Similar 
double logarithmic dependencies emerge in voltage- and 
temperature-dependent contributions to be analyzed be- 
low. 

Let us again stick to the limit eV ^ T . At voltages 
D/Ll < eF < D/Ll we obtain 



dl_ 

dV 



R 



£Th 

\eV\ 



1 



^Ly^^\eV\Ll 



,2 r2 



irgL 



D 



In- 



\eV\D 



(53) 

Note that in contrast to the expression for k (|52|) the 
voltage-dependent part of the interaction correction H53(l 
does not involve the upper frequency cutoff ~ l/r. At 
lower voltages eV < exh the differential conductance 
again saturates to the value (|49|) with k defined in Eq. 
(|52|l . The above results do not coincide with ones ob- 
tained in Ref. jl^ for the same physical situation. 

It is worth pointing out that even though the combi- 
nation Ly/gL in the above expressions always remains 
small, in some cases this smallness can be compensated 
by large logarithms, and the formal expansion in 1/g can 
become insufficient. In this situation the I — V curve can 
be evaluated with the aid of the following expression 



P(E) = 2^<5(i;)e-*(°) + 27r— e-*(^\ 

dE 



where 



■In- 



■In- 



Th 



(54) 



(55) 



This function also enters into the expression for the den- 
sity of states of normal 2d filmsiS 



(56) 



We also note that in the case of thicker films Lz ^ I our 
results are consistent with the expression for the density 
of states^. 



In the presence of massive gate electrodes close to our 
sample the fluctuating electromagnetic fields do not ex- 
tend outside its volume, and the term V* in Eq. Ipl?^ 
gets modified. In this limit we find 



dl 
dV 

where 



R 



eTh_ „ 
\eV\' 



^In^ 
ngL 



\eV\L 



D 



2Ly 

ngL 



In 



2 _S 
tD 



(57) 



(58) 



If one considers the limit T ^ eV, it suffices to sub- 
stitute 0.56r instead of eV in Eqs. (jST)) , and one 
will arrive at the expression for the linear Andreev con- 
ductance Ga{T) valid with the logarithmic accuracy. 



D. Effect of interactions: Quasi-ld structures 



Let us finally turn to quasi-ld structures L ^ Ly ^ Lz 
and examine the effect of electron-electron interactions 
on the Andreev conductance in the absence of screening 
gates. Similarly to the 2d case, there also exists the pa- 
rameter region Dq^ <C ui and aV q^ \og{l/ L^q^) ^ u; , 
which gives the leading contribution to J{t). We obtain 



00 00 



dq 



27r aTajq'^ 
|coth— [cos(ajt) — 1] — isin(cjt)| 
Since 1/q cannot exceed L, one should choose 



L-1 



(tF \n[crr/ujLl] 
Performing the ^-integration one finds 



(59) 



(60) 



J(t) 



duj 



1 + VujRC 



OJ 



! coth — [cos(wt) 



3in(cjt)| 



where 



C = 



L 



ln[i2/r] 



(61) 



(62) 



is the capacitance of a long cylinder. The expression 
(|61|l is rather similar to that obtained within the P{E)- 
theory—. There are at least two differences, however. 
One of them is that here the combination 1 4- \fujRC en- 
ters in the denominator of Eq. (|61|l . while in the modeU^ 



one has 1 



,2^2^2 



J' 



where Cj is the capacitance of a 
tunnel barrier at the NS interface. This capacitance is 
neglected in our analysis and can be trivially restored, if 



8 



^ 0.1 

< 



0.01 




FIG. 2: The zero temperature differential Andreev conduc- 
tance Ga = dl/dV of quasi-ld diffusive NS systems. We 
have set RCeTh = 0.01. The dimensionless parameter 8/g 
equals to 0, 0.1, 0.5, 1.1 (from top to bottom). 




|eV|/. 



FIG. 3: The (normalized) value of the interaction correction 

to the conductance SGa = Ga — G^2^ as a function of voltage 
for Ly = 0.3L, Lz — O.IL. One observes the saturation regime 
at eV ^ ETh, the Id logarithmic regime (EB at eV > 
which eventually crosses over to the 3d regime 1471 at eV ~ 
lO^eTh. 



needed. Another - more important - difference is that 
the frequency integral in Eq. H61|l cannot be extended 
below the Thouless energy. This low frequency cut-off 
results in the saturation of the interaction correction at 
small voltages and temperatures. This effect is not con- 
tained in the standard P(-E)-theory which does not keep 
track of quantum coherence of electrons entering an ef- 
fective environment and which also does not account for 
the sample geometry. 

In the limit T — > the asymptotic expression for the 
function J{t) at t » RC reads 

J{t) = - [In (eThi?C) - ci(|i|eTh) + isignt si(|t|eTh)] , 

. ^ CO . CO . 

where cit — — dz cos z/z, sit = — dz sin z/z. In 

particular, for t <^ l/exh we get 



Jit) 



In 



(64) 



where 7 ~ 0.577... is the Eulcr constant. At longer times 
t » l/exh the function J{t) approaches the constant, 
J{t) ~ — K, where 



Ki = — In 

g RCfTh 



8 e2p|,rin(LVr) 



= - In 
9 



(65) 



The I — V curve can now be easily evaluated. It is 
clear from Eq. H61|) that in the limit of high voltages 
eV ^ l/RC the response of an effective environment 
(normal metal) will be that of an i?C-transmission line. 
In this limit the interaction correction to the Andreev 
conductance will behave analogously to that studied for 
single electron tunneling within the P(i?)-theory^^. At 



lower voltages exh <C <C 1 /RC we find 
R 



d]_ 

dV 



-2Rl 

-87/3 



W\ 



(66) 



r(! + i) 



In the limit of large conductances 5 2> 1 this expression 
can also be rewritten in a simpler form 



dl 



R 
2if 



W\ 



1 



5 exh 



(67) 



At smaller voltages eV < exh we again observe the sat- 
uration of the differential conductance to the value l|49() 
where k is now determined by Eq. H65|) . 

The above results remain practically the same also in 
the presence of screening electrodes in the vicinity of the 
sample. In that case for g ^ 1 we again recover the 
logarithmic correction H67|l at large voltages while in the 
low voltage limit eV <C exh we find 



dl 
dV 



Re- 



32 f eV 



155 V^Th 



(68) 



As before, in order to obtain the linear Andreev conduc- 
tance Ga (T) in the regime T ^ eV with the logarithmic 
accuracy, in the above expressions it suffices to substitute 
0.56r instead of eV. 

The differential conductance of quasi-ld diffusive NS 
structures in the presence of electron-electron interac- 
tions is displayed in Fig. 2 for several values of the pa- 
rameter 8/g. At larger voltages one observes a competi- 
tion between two types of ZBA - conductance enhance- 
ment due to quantum interference of electronsiSiiiiiSii^ 
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and conductance suppression due to electron-electron in- 
teractions. For small conductances g < 16 the second ef- 
fect prevails and dl/dV increases with voltage. At large 
conductances g > 16 Coulomb effects become weaker 
and dl/dV shows the opposite trend. For smaller volt- 
ages eV < CTh the Andrcev conductance saturates to the 
voltage-independent value 149|) . (|65|l for all values of g. 
In this regime dl/dV is always suppressed by electron- 
electron interactions and this effect becomes stronger 
with decreasing g. 

Note, that the quasi-ld logarithmic behavior H67|) is 
rather robust and can manifest itself for a wide range of 
voltages even if the sample length L is not really much 
bigger than its width and thickness. An example of 
such a behavior is presented in Fig. 3 for the sample 
with Ly = 0.3L,Lz = O.IL. In this particular case the 
logarithmic dependence of Eq. H67|l remains applicable 
within a very wide voltage range exh ^ eV ^ lO^exh, 
while the 3d regime 1)47(1 may set in only at huge volt- 
ages eV > lO^exh- 



VI. CURRENT NOISE 

Our effective action analysis also allows to study fluc- 
tuations of Andreev current in diffusive NS structures in 
the presence of electron-electron interactions. Here we 
restrict our attention to the current noise. The noise 
spectrum is expressed in terms of the current operators 
/ in a standard way as 



S{t,t')^l{lit)Iit')+Iit')Iit) 



This expression can easily be translated into 



(69) 



6Sa 



(70) 

where SSa is defined in Eq. H27|l . Employing this gen- 
eral expression for the Andreev action together with the 
definitions of current and noise ((28(1 . ((70(1 one can derive 
the following relationship 

5(-,^,T) = eX:coth(^)/(|±y,r), 

.(71) 

which remains valid irrespective of the dimensionality at 
frequencies and voltages smaller than both the inverse 
charge relaxation time 1 /RC and the inverse elastic scat- 
tering time l/r. 

Note that the analogous expression between the cur- 
rent and the noise spectrum was derived in Ref. 
for normal tunnel barriers with interactions treated by 
means of the P(_E)-theory and in Ref. ,2^ within a more 
general theoretical framework. Here Eq. ((71(1 was derived 
for the Andreev current and the Andreev noise spectrum 
in diffusive NS structures in the presence of electron- 
electron interactions treated within a microscopic the- 
ory. Eq. I(71() has exactly the same form as thaiM^, one 



should only substitute 2e instead of e for the charge of 
the charge carriers. In the zero frequency limit Eq. 1(71(1 
reduces to the generalized Schottky relation 



S{uj,V,T) = 2ecoih{eV/T)I{V,T), 



(72) 



previously obtained in Ref. |25| in the absence of interac- 
tions. 

Combining Eq. ((71(1 with the results derived in the pre- 
vious section one can fully describe the effect of Coulomb 
interaction on Andreev current noise in diffusive NS 
structures. Hence, there is no need to go into further 
details here. Let us only mention that - similarly to 
the Andreev conductance - at voltages/frequencies be- 
low exh we expect shot noise to be suppressed by the 
factor exp(— k). At larger voltages eV > exh and in the 
low frequency limit the shot noise for quasi-ld NS struc- 
tures should scale as 



S cx \V\ 



8/3-1/2 



(73) 



VII. DISCUSSION 



The main observations of the present work can be sum- 
marized as follows. In hybrid NS structures electron- 
electron interactions generate two types of corrections 
to Andreev current. One of them is just a voltage- 
independent renormalization of the I — V curve by the 
factor exp(— k), while the other in general depends on 
voltage and/or temperature. The net effect of electron- 
electron interactions is always a reduction of the Andreev 
conductance. This reduction approaches its maximum 
value at temperatures and voltages below the Thoulcss 
energy of the sample. In this low energy regime the tun- 
nel barrier and a normal metal act as a single coher- 
ent scatterer, whereas at energies above exh they can be 
treated as independent ones in the spirit of the P{E)- 
theory. 

The magnitude of Coulomb effects essentially depends 
on the system dimensionality and is characterized by 
the dimensionless parameter k. In 3d systems we find 
K ^ l/p']?l'^ ^ 1, i.e. in this case Coulomb suppres- 
sion of the Andreev conductance is weak except for very 
disordered samples. In the 2d case n is proportional to 
the parameter l/pj^lLz (which is again small) but is en- 
hanced by the two logarithmic factors (cf. Eq. 1(52(1 '). 
i.e. the net effect of electron-electron interactions is not 
necessarily weak in this case. 

This effect becomes even more pronounced in quasi-ld 
systems in which case the value k is defined by Eq. 1(65(1 . 
Obviously, the Coulomb suppression of the Andreev con- 
ductance can be very strong for such systems. Combin- 
ing Eqs. (O and (gHI) with Eq. and having 
in mind that for typical values of the Fermi velocity in 
metals one has jvp ^ 1, we arrive at a very simple 
estimate for the Andreev conductance in the presence of 
electron-electron interactions: 



Ga/G^ ^ 1/N'^t 



(74) 
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where A^ch ~ Pp^ is the number of conducting channels 
in the normal wire. In Eq. H74|l for simplicity we have 
disregarded an additional weak (logarithmic) dependence 
on L and F, cf. Eq. (|^ . Within this accuracy, Eq. ifTHl 
demonstrates that Coulomb suppression of the Andreev 
conductance in quasi-ld diffusive samples is determined 
by the number of conducting channels to the power 8/g. 
We believe this prediction can be directly tested in future 
experiments. 

Recently the authorsiX have experimentally investi- 
gated the I — V curves for hybrid systems composed 
of multi-walled carbon nanotubes (MWNT) attached to 
a superconducting electrode. The measurements^^ per- 
formed below the superconducting critical temperature 
Tc have revealed power-law dependencies of the differ- 
ential conductance on voltage and temperature dl/dV oc 

and dl/dV cx F"^. At smaller voltages V < 14at 
1 mV the conductance was found to saturate to a con- 
stant value. Also at T > Tc the power-law dependencies 
have been observed dl/dV oc T"" with the value aN 
somewhat smaller than as- 

It is interesting to compare our theoretical predictions 
with these experimental results. The effect of conduc- 
tance saturation observed in Ref. just matches with 
our prediction l|7I|l and the value eT4at appears to be in a 
good agreement with the Thouless energy exh estimated 
for MWNTii. Also the observed power-law dependen- 
cies are consistent with our theoretical picture. Identify- 
ing as = 8/g — 1/2 and making use of the experimental 
value as — 1.28, we can estimate the effective dimension- 



less conductance for MWNT as g « 4.5. Above Tc for the 
superconducting electrode one also expects the power-law 
dependence of the conductance on temperatur o^^i^^ . Ne- 
glecting the impedance of the bulk electrode one would 
trivially get aN = 2/g « 0.45 which is slightly smaller 
than the measured values ajv ~ 0.77 close to Tc and 
ajv ~ 0.55 at higher temperatures. It is clear, how- 
ever, that above Tc the impedance of the superconduct- 
ing electrode Rs becomes non-zero and should also be 
taken into account. Accordingly, the value a^v should be 
modified as aAr = 2/^ + 2/gs, where gs = Rq/Rs- This 
additional contribution to aN could further improve the 
quantitative agreement between our theory and the ex- 
perimental observationsii. 

Finally, we would like to point out that it would also be 
interesting to experimentally study the effect of Coulomb 
interactions on the shot noise in quasi-ld hybrid NS sys- 
tems, like, e.g., ones investigated in Ref. ^3 Accord- 
ing to the results obtained in Sec. VI, at low voltages 
one expects suppression of the shot noise by the factor 
^ ^cii ^tiils higher voltages the power-law depen- 
dence (|75|l should be recovered. Simultaneous current 
and shot noise measurements could help to experimen- 
tally verify the general relation (|71|l between the noise 
spectrum and Andreev current in hybrid NS structures. 
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